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Abstract 

We face the problem to determine the slope dependent current during the epitaxial growth process of a crystal 
surface. This current is proportional to 5 — p+ — p_, where p± are the probabilities for an atom landing on a 
terrace to attach to the ascending (p+) or descending (p_) step. If the landing probability is spatially uniform, the 
current is proved to be proportional to the average (signed) distance traveled by an adatom before incorporation 
in the growing surface. The phenomenon of slope selection is determined by the vanishing of the asymmetry 8. We 
apply our results to the case of atoms feeling step edge barriers and downward funnelling, or step edge barriers 
and steering. In the general case, it is not correct to consider the slope dependent current j as a sum of separate 
contributions due to different mechanisms. 

Resume 

Courants de surface et selection de la pente en croissance cristalline La croissance epitaxiale d'une 
surface cristalline peut etre caracterisee par un courant de surface J, dont la partie j qui depend de la pente est 
etudiee. Celle-ci est proportionnelle a S — p+ —p~, ou p± sont les probabilites qu'un atome depose sur une terrasse 
se colle a la marche montante (p+) ou descendante (p— ). Si la probability spatiale d'atterrissage est uniforme, le 
courant est aussi proportionnel a la distance moyenne (avec signe) parcourue par chaque atome. Le phenomene 
de la selection de la pente est determine par la condition 5 = 0. Les resultats ainsi obtenus sont appliques aux cas 
barriere de marche plus downward funnelling et barriere de marche plus braquage (steering). Dans le cas general, 
le courant j ne peut pas etre considere comme la somme de contributions separes dues aux differents mecanismes. 
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1. Introduction 



Molecular Beam Epitaxy is a well known and widespread technique to growth layers of metal and semiconductor 
crystals. The growth process of a high symmetry crystal surface can be described through a ballistic flux of 
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particles impinging on the growing surface and the following thermal diffusion of adatoms, which finally attach 
to a preexistent step or nucleate new terraces [1,2,3]. 

One crucial aspect of the growth process is its possible unstable character, due to deterministic mechanisms 
which prevent the growing surface to remain flat. In this paper we focus on mound formation and we will discuss 
its description via a mesoscopic surface current J, which enters in the evolution equation dtz = Ft — d x J, where F 
is the flux of incoming particles and z is the local height. Experiments show that dynamics can produce coarsening, 
with the formation of mound facets of constant slope m* [4] . 

In the theoretical seminal papers [5,6] the cause for mound formation was identified in the existence of step edge 
barriers, i.e. an additional barrier hindering the descent of steps. Since then, many efforts [2] have been devoted to 
the formalization of this idea and to the derivation of quantitative predictions for the evolution of mounds. These 
efforts have often combined phenomenological approaches with attempts of a rigorous derivation of J starting 
from the microscopic dynamics of adatons. 

The main piece of the current J is the so called slope dependent current j(m), depending on the local slope 
m of the surface. J also contains terms depending on higher order derivates, but the rising of the instability and 
the possible formation of mounds with a constant slope m* only depend on j(m). Therefore, it is natural that 
special attention has been devoted to its determination. Analytical approaches have followed two mainstreams: 
coarse-graining procedures [7,8,9] to pass from step-dynamics to mesoscopic dynamics (see next Section) and the 
evaluation of j(m) through the average (signed) distance [10,11] walked by adatoms before being incorporated in 
the growing crystal (see Section 5). 

A recent paper by Li and Evans [12] has renewed the interest on the slope dependent current. Authors claimed 
that standard phenomenological continuum theories are inappropriate to describe mound slope selection. After- 
wards some of their claims have been corrected [13], but their work has shown that j(m) should be evaluated 
with great care. For these reasons, in the following we reconsider the problem, giving a formulation as general 
as possible for the slope dependent current and for the evaluation of the selected slope m* , determined by the 
condition j(m*) = 0. We prove that the current j(m) is proportional to the asymmetry 5(L) = p+(L) — P-(L) 
between the probabilities for an atom deposited on a terrace of size L to be incorporated into the upper (p+) 
and lower (p_) steps. This proof does not refer to any specific microscopic mechanism. An important by-product 
is that, in general, j(m) can not be considered as a sum of contributions due to separate mechanisms. Finally, 
we also discuss the alternative formulation for the current in terms of the averaged (signed) distance walked by 
adatoms before incorporation. 



2. The current 



In the case of conserved growth (no desorption, no overhangs), a useful concept to study the dynamics of the 
surface is the mesoscopic current J, entering via the evolution equation dtz = Ft — d x J for the local height z. It 
is worth stressing that the dynamics of the surface, and therefore the current J, are determined by the dynamics 
of steps: adatoms enter only through their attachment (and detachment) rate to steps. 

A piece of surface of slope m looks differently, according to the value of the slope. For large m (Fig. la) we have 
a sequence of all uphill or downhill terraces of size L ~ l/|m|. For small m (Fig. lb) we have a mix of different 
types of terraces of size L ~ £d, where to is the nucleation length [2]: in this case, the average slope is determined 
by uncompensated uphill and downhill terraces. The two pieces of surface look differently because nucleation on 
terraces prevent them being larger than Id. 

The analytical determination of the current J starts from the large slope case (see Fig. la). In a one dimensional 
picture, the flux FL of atoms landing on the terrace splits in two currents j±(L) = FLp±(L), where p±(L) are 
the probabilities that an atom attaches to the ascending (p+) or descending (p~) step. The velocity of the n—th 
step is simply equal to \j-(L n +i) + j+(L n )]. The sum (p+ +p~) = 1, while 8 = (p+ —p~) defines the asymmetry 
S. 

Following the method introduced in Ref. [7], the displacement of step n during the deposition of one monolayer 
can be approximated as 

x n (t + l/F) - x n (t) = ~\j-{L") +3+(L% (1) 
where L' = i(L„ + L„_i) and L" = ±(L„ + L n +i). It is useful to sum the quantity L n to both sides, so as to have 
x„{t + 1/F) - x n (t) + L n = -\\L" + L' - 2L n + L'S(L') - L"5(L")]. (2) 
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Figure 1. (a,b) Schematic of a growing one dimensional surface at (a) large and (b) small slope (figure b has been shrinked). 
(c,d) Coordinates for the models studied in Section 3: (c) Downward Funnelling and Ehrlich-Schwocbel effect; (d) steering and 
Ehrlich-Schwocbcl effect. 



The quantity on the left can be approximated as { — l/mF)dtZ = {l/mF)d x J, where m 
quantity on the right can be worked out using the relations [7,14] 



L" +L' - 2L„ 



[L d xx L + L(d x Lf 



L'S(L') - L"S(L") « [l! - L")d L [L5{L)\ « Ld x Ld L [LS(L)]. 
Finally, we get 



-^-=d x J = —d x 
rat m 



Ld x L LS(L) 



4 2 

so that the total current J comes out to be 



J = \FL5{L) - \FLd x L. 



1/L is the slope. The 

(3) 
(4) 

(5) 

(6) 



If we use the slope m, the second term takes the form jF(d x m)/m 3 , which was already found in Refs. [7,8]. In 
the following we will focus on the slope dependent part, 



j = \FLS{L). 



(7) 



We stress that Eq. (7) is the most general form, not depending on any assumption on the miscroscopic processes 
occurring at the surface. It is valid for large slope, m > 1/£d (Fig- la). In the opposite limit, m < 1/£d (Fig. lb), 
j(m) is linear [7] in m, as also expected from symmetry reasons for a slope dependent current: j(m) = j{l/lD)t-Dm, 
where j(l/£o) = ^F£d5(£d) is evaluated according to (7). 

The condition of instability of the flat surface reads j'(O) > 0, i.e. j(l/£n) > 0, as trivially shown by a linear 
stability analysis [15]. Finally, if mounds develop facets with constant slope m* , the current must vanish on it: 
j(m*)=0. 



3. Downward funnelling, step edge barriers and steering 



Let us now consider the following model (see Fig. lc). Atoms deposited within a distance c from the descending 
step are incorporated in it (downward funnelling, DF), while atoms deposited in the remaining (L — c) portion 
of the terrace diffuse freely, feeling an additional (Ehrlich-Schwoebcl, ES) barrier at the descending step. Atoms 
deposited in the c region give a trivial contribution Fc to j-. As for the others, we must solve the diffusion 
equation F + Dp"(x) = for c < x < L and p"(x) = for < x < c, with boundary conditions p'(0) = p(0)/£es, 
p(L) — 0, and p(x),p'(x) continuous in x = c. The quantity £es > is the well known Ehrlich-Schwoebel 
length and measures the additional barrier felt by an adatom in the sticking process to the descending step. It is 
straightforward to get 

> m _ F (L-cf , F (L-c)(L + c + 2£ ES ) 
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Figure 2. (a) The selected slope cm* as a function of A, Eqs. (12-13). (b) Full line: the selected slope cm* as a function of Ies/c 
(DF+ES model, Eq. (14)). Dashed line: the slope m* such that j(m*) — for the stecring+ES model (Section 3). There is no slope 
selection in this case, because j'(m*) > 0. 

The contributions of atoms deposited in the (L — c) region to the currents j± are D\p'(L)\ and Dp'(0), respec- 
tively. Therefore 

S -=Fc + ¥£f = FI* A L), i+ = f C£^ + £±^ sF ^ (£) . (9) 
2 L + Ies 2 L + «ES 

Their sum gives the total flux of particles arriving on the terrace: j+ + j- = FL, while their (semi-)difference 
gives the current j = ^FLS(L) = \ — j-), i.e. 

. _ F l E s(L - 2c) - c 2 _ F £ E s(l - 2mc) - mc 2 

3 ~ 2 L + Ies ~ 2 1 + mt E s ' ( ' 

If c = 0, we find the well known result j = \FLZ E s/{L + £es) = \FZ E s/(l + m£ E s)- 

Eq. (10) shows that it is not generally possible to write j as a sum of two separate contributions, j = jES+jDF, 
due to the Ehrlich-Schwoebel (ES) effect and to the downward funnelling (DF), respectively. In fact, it is instructive 
to consider the two limits Ies = and Ies = oo. For £es = oo, = Fc and j + = F(L — c), i.e. j- is entirely 
due to downward funnelling and j + to the ES effect. It seems natural to write j = j DF + Jes, with j DF — —\Fc 
and jes = \F(L — c). For £es = we are induced to assume j = jdf = —Fc 2 /2L. Therefore, if we want to write 
j = jDF +jES, the expression for jdf depends on £es- it is constant (not depending on m) for large barriers and 
it is linear in m for weak barriers. The conclusion is that j(L) should be handled as a whole. 

Let us now consider a situation where the flux of particles on the terrace is not uniform, see Fig. Id. Because 
of steering effects [16], we assume that all atoms destined to land within distance c from the ascending step are 
steered and finally land in the c region close to the descending step (of the upper terrace) . So, this region undergoes 
an effective flux 2F. Working out calculations [17] similar to those made here above, we get the result 

. = F £ ES L - 2c(L -c) = F 2c 2 m + £ ES - 2c 

3 2 L + £ ES ~ 2 l + m£ E s ' ( ' 

This current is always positive if £es > 2c, or it changes sign in m* = (1 — £es /2c) / 'c, if £es < 2c (see Fig. 2b, 
dashed line). However, in the latter case j'(m*) > 0, which implies instability of the slope m* [18]. So, in this 
model there is no slope selection: either the surface is always unstable or it is metastable [19] . 



4. Slope selection 

First, let us apply the condition j(m) = to a model studied by Li and Evans [12]: 

j+=F(L~c)P+, j-=Fc + F(L-c)P- (12) 
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If P± are the probabilities for atoms deposited in the (L — c) region to attach to the step, this model is similar 
to what we studied in the previous Section. For constant P±, the condition j + — j_ , i.e. 5 = 0, gives the following 
expression for the selected slope m* , 



m * = <TTA) ' (13) 

where A = P+ — P_. Note that A is not the total asymmetry 5, which vanishes for m = m* , but the asymmetry 
for atoms deposited in the (L — c) portion of the terrace only. In Fig. 2a we report cm* (A). Eq. (13) perfectly 
matches the numerical results given in [12]. 

Let us now turn back to Eq. (10). In this case, we get 

(14) 



c(c+2£ ES )' 

In Fig. 2b (full line) we plot cm* as a function of £es/c. For small and large £es/c, m* = £es/c and m* = l/(2c), 
respectively. Therefore, for large ES barrier the selected slope corresponds to a terrace size equal to 2c, as expected 
by a trivial compensation of DF and ES effects. 

The use of (10), as well as of (12), to find m* is limited by the constraint m* > 1/£d- For Eq. (10) this means 
£es > c 2 /(£d — 2c) ~ c 2 /£r>; in the opposite limit, £es < (? /£d, 3(X/£d) < and no selected slope exists: this 
happens because the DF effect is so strong to induce a stabilizing (downhill) current at all slopes. 



5. Another expression for the current 



The expression j(L) = ^FLS(L) is the most general one for the slope dependent current on a region of (large) 
slope m — 1/L. The quantity S(L) = (p + — p_) measures the asymmetry between the sticking probabilities to the 
upper and lower steps. The flux F is usually assumed to be spatially uniform, apart from fluctuations. However, as 
anticipated in the previous Section, because of atom-substrate interactions, steering effects may occur and atoms 
are no more deposited uniformly on the terrace. In spite of this, the espression j(L) — ^FLS(L) still continues to 
be correct. In the Introduction we mentioned a possible different evaluation of j(m), through the average signed 
distance d walked by an adatom, from the deposition to the incorporation site. In the following we are proving 
the equivalence of the two expressions, if the flux is uniform. If it is not uniform, the average distance d is not an 
appropriate quantity to determine j. 

The convention is to take d positive if the atom attaches to the ascending step (see Fig. 1): 

L L L 

d = ^ J dx[(L — x)p+(x) — xp-(x)] = J dx p+(x) — i J dx(p+(x) + p-(x))x, (15) 



with p±(x) being the probabilities that an adatom deposited in x attaches in x — (p-(x)) and in x = L (p+(x)). 
Since p± — (1/L) dxp±(x) and [p+{x) +p-(x)) = 1, we get d = Lp+ — = ^L8(L), so that 

j(L) = \FL8{L) = Fd. (16) 

In a Kinetic Monte Carlo simulation, we can easily implement the above expression and write 

JKMC = F — (17) 

where N r j are the total hops in the uphill and downhill direction, and N a is the total number of deposited atoms. 
Formula (17) was firstly introduced in Ref. [10] and has the merit to be valid for a surface of any slope. At author's 
knowledge, a rigorous derivation and a comparison with the mesoscopic current were missing. 

It is worth stressing that Eqs. (16-17) are no more applicable if the landing probability is not spatially uniform. 
In this case, j(L) is no more proportional to d, as shown by a trivial example: no downward funnelling (c = 0) 
and infinite ES barrier, £es = oo. The resulting current j = \FL can be written as j = Fd if the flux is uniform: 
in that case, d = \L. Different spatial distributions of the landing atoms modify d and therefore the expression 
j — Fd, but do not modify the correct expression j — hFL5, because the step dynamics would be unchaged. 
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6. Final remarks 



In this paper we mainly focused on the slope dependent current and the slope selection mechanism. We have 
shown that a correct derivation of this current is possible and its general expression has been found. The condition 
j(m*) = determines the selected slope. More generally, the shape z s (x) (or m s (x) = z' a (x)) of stationary solutions 
depends on the vanishing of the full current J in all the points, J(m s (x)) = 0. It is worthnoting that the conditions 
j(m*) = and J(m s (x)) = 0, when applied to a discrete or discrete-continuous model, are valid only averaging J 
on time scales not smaller than 1/F [12]. 

With regard to additional terms in the current, our derivation in Section 2 shows that a symmetry breaking 
term, having the form J = d x A(m 2 ), appears naturally when we coarsen step dynamics. This term is the only 
term surviving in a plain vicinal surface without any additional microscopic mechanism (no step edge barriers, no 
nucleation, no thermal detachment, no downward funnelling): so, in this sense, it should be considered the most 
fundamental one. Finally, J should also contain at least a Mullins-like term, J ~ d xx m, which may be due to 
several mechanisms [20]. 

However, even if the rigour in the derivation of the mesoscopic current has improved in the course of time, 
the evolution equation dtz = Ft — V • j is not much more than a phenomenological equation, specially in two 
dimensions where additional problems linked to step edge diffusion exist. The dynamics of a truly vicinal surface, 
which can be studied with much more rigour [21], shows that the full nonlinear equations governing the growth 
process of a real system are far more complicated. 
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